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Rotation of the NGC 1333 IRAS 4A2 Protostellar Jet 
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ABSTRACT 

The bipolar jet of the NGC 1333 IRAS 4A2 protostar shows a velocity gradient in the direction perpendicular to 
the jet axis. This lateral velocity gradient can be seen throughout the jet imaged in a silicon monoxide line, 2500-8700 
AU from the driving source, and is consistent with the rotation of the accretion disk. If this gradient is caused by the 
rotation of the jet around its axis, the average specific angular momentum is about 1.5 x 10^^ cm^ s^^. Comparison of 
the kinematics between the jet and the disk suggests that the jet-launching region on the disk has a radius of about 2 AU, 
which supports the disk-wind models. The angular momentum transported away by the jet seems to be large enough 
for the protostar to accrete matter from the disk, confirming the crucial role of jets in the early phase of star formation 
process. 



Subject headings: ISM: individual objects (NGC 1333 IRAS 4A2) 
stars: formation 



ISM: jets and outflows — ISM: structure 



1. INTRODUCTION 

In star formation processes, from a dense core of interstel- 
lar cloud to a main-sequence star, the specific angular momen- 
tum has to decrease million-fold (Bodenheimer 1995), and the 
loss of angular momentum is essential. The spin-down mech- 
anism is not well understood, but the transportation of angular 
momentum by a protostellar jet is thought to be one of the im- 
portant processes (Ray et al. 2007). The jets are thought to be 
driven by rotating disks through magneto-centrifugal processes 
(Blandford & Payne 1982; Ray et al. 2007; Pudritz et al. 2007). 
However, quantitative understanding of the jet-driving mecha- 
nism has been difficult because examples showing rotation in 
both disk and jet are rare. To understand this process quantita- 
tively, it is important to study the kinematics of a system show- 
ing the rotation in both disk and jet. 

NGC 1333 IRAS 4A2 is a Class protostar that belongs to 
a binary system in the Perseus star-forming region at a distance 
of 235 pc from the Sun (Sandell et al. 1991; Lay et al. 1995; 
Looney et al. 2000; Hirota et al. 2008; Enoch et al. 2009). IRAS 
4A2 drives a prominent jet seen clearly in many molecular lines 
(Blake et al. 1995; Lefloch et al. 1998). The IRAS 4A2 bipolar 
jet imaged in an SiO line is extremely well coUimated and suit- 
able for a detailed study (Choi 2005). Comparison of the proper 
motion from molecular hydrogen images with the line-of-sight 
velocity from the SiO spectra suggests that the jet axis is close 
to the plane of the sky with an inclination angle of ~10?7 (Choi 
2005; Choi et al. 2006). 

It was previously noticed that the SiO jet displays a lateral 
velocity gradient, which was puzzling (Choi 2005). While rota- 
tion is an obvious explanation, detailed interpretations were de- 
ferred until the rotation direction of the circumstellar disk can 
be determined observationally. The disk of IRAS 4A2 is bright 
in the NH3 lines and elongated in the direction perpendicular to 
the jet axis (Choi et al. 2007). Recent observations showed that 
the disk displays a Kerplerian-like rotation (Choi et al. 2010). 
The derived mass of the central protostar is 0.08 Mq, and the 
collapse age is ^50,000 yr The sense of disk rotation is con- 
sistent with the velocity gradient of the SiO jet, which strongly 
suggests that the jet is rotating around the flow axis. 

In this Letter, we present the results of our data analysis of 
the NGC 1333 IRAS 4A2 SiO jet. We describe the data and 



results in Section 2. In Section 3, we discuss the rotation of the 
IRAS 4A2 jet and its implications. 

2. DATA AND RESULTS 

Details of the observations and the results were presented by 
Choi (2005). The NGC 1333 IRAS 4 region was observed using 
the Very Large AiTay (VLA) in the SiO u = 0J = 1^0 line. 
The resulting image has a restoring beam of FWHM = 1'.'96 
with a natural weighting. When a higher resolution is needed 
for the analysis, an image made with a robust weighting was 
used, which has a beam size of 1'.'50. The natural weighting 
gives a higher signal-to-noise ratio, while the robust weighting 
gives a smaller beam size. 

There are several SiO emission peaks distributed along the 
bipolar jet of IRAS 4A2 (see Figure 2(a) of Choi 2005). The 
downstream half of the northeastern jet (SiO outflow peaks 1- 
3) is strongly disturbed by the interaction with an ambient cloud 
core (Choi 2005; Baek et al. 2009). In this Letter, we concen- 
trate our attention to the undisturbed part of the jet (Figure 1; 
SiO outflow peaks 4-12). 

To investigate the kinematics along the lateral direction (i.e., 
perpendicular to the jet axis) position-velocity (PV) diagrams 
were made along several cuts passing through the emission 
peaks (Figure 2). The reference position of each cut shown in 
Figure 1 is the peak position of the map integrated over the 
whole velocity interval. The exact location of the reference po- 
sition, however, is not important because all the calculations 
in the discussion require position differences, not the absolute 
position. All the PV diagrams show that the SiO emission is 
distributed along diagonal lines. The sense of velocity gradi- 
ent is consistent throughout the jet: the northwestern side is 
blueshifted relative to the southeastern side. 

3. DISCUSSION 

3.1. Jet Rotation 

Since the sense of lateral velocity gradient is consistent 
throughout the whole system (the southwestern jet, the disk, and 
the northeastern jet), the most natural and obvious explanation 
is that the velocity gradient is caused by the rotation around 
the jet/disk axis. Alternative explanations for the lateral veloc- 
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[See http://minho.kasi.re.kr/Publications.html for the original high-quality figure.] 

Fig. 1. — Color composite image of tlie NGC 1333 IRAS 4A2 bipolar 
jet in the SiO u = J = 1 ^> line (see Figures 2 and 3 of Choi 2005). 
Blue, cyan, green, yellow, and red images show the emission at Vlsr = 
-9.4, -4.7, 0.0, 14.8, and 19.6 km s"^, respectively, with a bandwidth 
of 4.7 km s~^. Shown in the bottom left-hand corner is the restoring 
beam: FWHM = l'.'96. Straight lines: cuts for PV diagrams (Figure 2). 
The PV cuts are perpendicular to the jet axis. A short tick marks the ref- 
erence position on each cut. Solid dots: the 3.6 cm continuum sources 
(Reipurth et al. 2002). 



ity gradients of outflow were proposed previously (e.g., Soker 
2005; Cerqueira et al. 2006), but they do not seem to apply to the 
case of IRAS 4A2 (see the discussion in Section 4.3 of Chrysos- 
tomou et al. 2008). The linear pattern of emission in the PV di- 
agrams (Figure 2) suggests that either the SiO emission comes 
from a cylindrical layer around the flow axis or the SiO jet is 
rotating like a rigid body. At several positions (e.g., cuts 10 and 
11) the SiO emission is stronger at both ends of the diagonal 
than at the middle part, probably owing to the limb-brightening 
effect. It should be noted that some emission peaks are off the 
diagonal lines (e.g., the -8.1 km s^^ peak of cut 7), suggest- 
ing that the jet structure is more complicated than our simple 
description. There is also some ambiguity in selecting the end 
points of the diagonal lines at a few cuts (e.g., the red ends of 
cut 7 and cut 10). The uncertainty introduced by this ambiguity 
does not significantly affect the conclusions in Sections 3.2 and 
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[See http://minho.kasi.re.lcr/Publications.html for the original high-quality figure.] 

Fig. 2. — PV diagrams of the SiO line along the cuts perpendicular to 
the jet axis (Figure 1). The horizontal axis is the angular distance from 
the reference position shown in Figure 1 . The PV diagrams of cuts 4, 
5, 8, 10, 11, and 12 are from a map made with a natural weighting, 
and the rms noise is 0.22 K. The PV diagrams of cuts 7 and 9 are 
from a map made with a robust weighting, and the rms noise is 0.40 K. 
Contour levels start from "io and increase by la, where o is the rms 
noise. Dashed contours are for negative levels. Dotted straight lines: 
position-velocity gradient. 



3.3. 

To understand the jet rotation quantitatively, several physical 
parameters are plotted as a function of distance from the driv- 
ing source (Figure 3). The jet radius, Rj, is the half-distance 
between the end points of the diagonal lines shown in Figure 2. 
The rotation speed, v^, is half of the velocity difference between 
the end points of the diagonal lines, corrected for the inclination 
angle. The angular speed is O = v^/Rj, and the specific angu- 
lar momentum is L = RjV^. These parameters were calculated 
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Fig. 3 — Properties of the IRAS 4A2 SiO jet. The horizontal axis 
shows the angular distance from IRAS 4A2. (a) Radius of the jet at 
each lateral cut. Names of the cuts are labeled. Cuts 4 and 5 (open cir- 
cles) are on the northeastern jet, and cuts 7-12 (filled circles) are on 
the southwestern jet. (b) Rotation speed, (c) Angular speed, (d) Spe- 
cific angular momentum. Dotted lines: linear fits to Rj, v^, and L. For 
n, straight lines do not give a good description of the data. Dashed 
curve: ratio of the linear fits to V4, and Rj . 
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Fig. 4. — Simple schematic of the IRAS 4A2 disk-jet system. The hor- 
izontal axis shows the distance along the jet axis from the protostar at 
the center of the accretion disk. The SiO-emitting layer of the jet is a 
cone with an opening half-angle of ~3?1. Loci of particles in the jet 
are shown in solid curves on the near surface and in dotted curves on 
the far surface. The pitch angle is ~2?9 in the detected part of the SiO 
jet. The jet may rotate faster in the upstream part where the SiO emis- 
sion was not detected. The jet is bipolar, one ejected to the right and 
the other to the left, but only one is shown. The counter-jet would be 
a mirror image of the jet shown here. The disk has a radius of 3 10 AU 
(Choi et al. 2010). The outflow foot-ring of 2 AU is drawn around the 
center of the disk, but it is very small in this scale and appears like a 
point. 

jet (Figure 3(b)). The angular speed also decreases with the dis- 
tance: the slope is relatively steep at small distances and then 
becomes shallow at large distances (Figure 3(c)). These behav- 
iors of U0 and fl are expected from the conservation of angular 
momentum when the jet has a finite opening angle. However, L 
increases with the distance (Figure 3(d)), which implies that ei- 
ther the angular momentum injected into the jet may vary with 
time or the correction factor varies as the jet propagates (see the 
discussion in Section 3.2). In any case, the fact that L is not 
decreasing suggests that the loss of angular momentum through 
the interaction with the ambient medium is not an efficient pro- 
cess. 

It is interesting to note that L of the paired data points at 
the same distance from the driving source (cuts 5/7 and cuts 
4/8) agree within the uncertainties. This physical connection 
between the opposite lobes of the bipolar jet may be natural 
because the jet materials at each pair were ejected at almost the 
same time from the jet engine (disk). Therefore, the physical pa- 
rameters of a specific part of the jet can be used to understand 
the properties of the jet engine at the time of ejection. 



assuming that the SiO emission at each lateral cut (cross sec- 
tion perpendicular to the jet axis) comes from a hollow ring. 
Otherwise, a correction factor may be applied to each parame- 
ter. For example, if the SiO emission comes from a uniform- 
density disk with a rigid-body rotation, the average rotation 
speed would be shown in Figure 3(b) multiplied by a factor 
of 2/3, and the average specific angular momentum would need 
a correction factor of 1/2. A more complete description of these 
parameters would require information on the jet structure such 
as rotation curve, density profile, and molecular abundance pro- 
file as functions of radius from the jet axis, but these details are 
unavailable and beyond the scope of this Letter. Multi-species 
multi-transition observations in the future can be helpful. 

Figure 3(a) shows that Rj increases almost linearly, and ex- 
trapolating the best-fit line to the driving source (position offset 
of zero) gives a small radius (100 ± 70 AU), which suggests 
that the shape of the SiO jet is similar to a cone with a finite 
opening angle. Figure 4 shows a schematic description of the 
jet structure. The rotation speed seems to decrease with the dis- 
tance, but the scatter is large in the downstream part of the SiO 



3.2. Jet Engine 

One of the interesting quantities in many models of jet driv- 
ing mechanism is the location of the jet-launching region on 
the disk. Simply extrapolating Rj to a small distance would not 
give the answer because the intrinsic shape of the jet stream 
line is unknown. Instead, the conservation laws of magneto- 
hydrodynamic winds can be used (Pudritz et al. 2007). The en- 
ergy and angular momentum extracted from the disk are in an 
electromagnetic form at the base (or foot) of the jet and almost 
completely converted to a kinetic form in the observed part of 
the jet (Anderson et al. 2003). Therefore, Rj and measured 
from the SiO jet can be used to derive the angular speed at the 
base of the jet, Oo> which can be converted to the radius of out- 
flow foot-ring, Rf. 

The IRAS 4A2 jet was seen in a molecular hydrogen line, 
and its proper motion was measured (Choi et al. 2006), which 
gives the velocity of the jet in the plane of the sky to be 
about 71 km s~^. If this value can be applied to the SiO jet, 
the total specific energy and angular momentum can be calcu- 
lated at each position. These quantities are conserved along the 
field line as protostellar outflows are considered cold magneto- 
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hydrodynamic winds (Mestel 1968; Pudritz et al. 2007). Then 
can be obtained by solving equation (4) of Anderson et al. 
(2003), which can be converted to the foot-ring radius of the jet, 



(1) 



where M* is the mass of the central protostar. Figure 5(a) shows 
Rf at the time of ejection for each lateral cut on the jet, which 
is distributed in the range of 1.5-2.8 AU. The average foot-ring 
radius is 2.0 AU. Figure 5(a) also shows that Rj is slowly in- 
creasing with the distance or jet propagation time. If this trend 
is real, it could mean that the outflow launching region on the 
disk is shrinking in size. 

The jet, however, is not necessarily launched from a thin 
ring, and the distribution of emission needs to be considered. 
For example, the emission of cut 10 is mostly concentrated at 
the two ends of the diagonal line in the PV diagram (Figure 2), 
and the Rf given in Figure 5(a) is a good representative value. 
By contrast, for cut 9, strong emission comes from the mid- 
dle part of the diagonal (Figure 2). If we take the two strongest 
peaks of cut 9 and repeat the same calculations, R j can be as 
small as 0.6 AU. Another implicit assumption is the Keplerian 
disk. As mentioned by Choi et al. (2010), the assumption of 
Keplerian rotation of accretion disk needs to be tested with fu- 
ture observations. Therefore, the uncertainties of i?/ and related 
quantities may be larger than the statistical uncertainties shown 
in Figure 5. 

The measured foot-ring radius of ~2 AU seems to favor the 
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Fig. 5. — Properties of the IRAS 4A2 jet engine. The horizontal axis 
shows the time of jet propagation from the driving source, assuming a 
constant proper motion of 0'.'064 yr"'^ (Choi et al. 2006). (a) Foot-ring 
radius inferred from the SiO jet properties at each lateral cut. Dotted 
Une: Unear fit. The slope is (1.9 ± 1.0) x 10"^ AUyr"\(b) Mass ejec- 
tion efficiency inferred from {Rj/Ra)^- The (Icr) uncertainties shown 
here include statistical uncertainties only. The uncertainty of the proto- 
stellar mass (0.08 ± 0.02 Mq) affects these values systematically, at a 
level of 10% for Rf and at a level of 20% for fm- Note that the quanti- 
ties shown here are calculated from the end points of the diagonal lines 
in the PV diagrams (Figure 2), which is suitable when the emission is 
concentrated at the end points. Otherwise, a range of values should be 
considered. For example, if we take the two strong peaks of cut 9 near 
the middle part of the diagonal line (Figure 2), Rf can be as small as 
0.6 AU, and fm can be as high as 4%. 



disk- wind models (Pudritz et al. 2007) that predict jets launched 
from an extensive area on the disk with a radius of a few AU. 
However, jets coming from the inner edge of the disk at a few 
stellar radii (Shu et al. 2000), cannot be ruled out, because it is 
possible that there is yet another component of jet that is not 
traced by the SiO line. Future observations of the jet rotation 
with atomic tracer lines would provide a more complete answer. 

3.3. Angular Momentum Transport 

If enough mass can be loaded, the jet can carry away the 
angular momentum in the disk so that the mass accretion can 
continue. The required ratio of mass ejection to mass accretion 
can be approximately obtained by comparing the specific an- 
gular momenta in the jet and disk (Pelletier & Pudritz 1992; 
Anderson et al. 2003; Pudritz et al. 2007). The mass ejection 
efficiency can be calculated by 



_M^^fRl 



(2) 



where is the mass outflow rate. Ma is the mass accretion 
rate, and Ra = (L/fio)^^^ is the Alfven radius. For the SiO jet 
of IRAS 4A2, fm is only a few percent (Figure 5(b)), which is 
slightly smaller than the typical value (5-10%) of young stellar 
objects (Ray et al. 2007). Therefore, the jet seems to be playing 
an essential role in the growth of the protostar IRAS 4A2. 

While rotating jets of some T Tauri stars were observed pre- 
viously (Bacciotti et al. 2002; Ray et al. 2007; Launhardt et 
al. 2009), they are young stellar objects well past the main ac- 
cretion phase, and the transport of angular momentum is less 
critical. Rotation of the HH 26 flow driven by a Class I source 
was reported by Chrysostomou et al. (2008), and the derived 
foot-ring radius is 2-4 AU. (Also see Zapata et al. (2010) for 
a possible rotating jet of a young object in an ambiguous evo- 
lutionary stage.) It has been elusive to detect a rotating jet of a 
protostar in the very early (Class 0) phase of evolution when the 
removal of angular momentum is crucial. For example, detec- 
tions of lateral velocity gradients were reported for the HH 211 
jet (Lee et al. 2007), but it was later found that the sense of gra- 
dient can be reversed at other parts of the jet (Lee et al. 2009), 
which suggests that the lateral velocity gradients in this jet are 
not caused by rotation. By contrast, the IRAS 4A2 Class sys- 
tem presented in this Letter shows the velocity gradient very 
consistently throughout the bipolar jet, which provides reliable 
measures of the rotation kinematics and allows a test of models 
of jet driving mechanism. 
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